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ABSTRACT 

We report on a 50 ksec XMM-Newton observation of the double pulsar system J0737-3039 performed 
on April 2004. We present results of the spectral analysis of these data combined with the much 
shorter Chandra pointing performed on January 2004. Black body emission with effective temperature 
of 0.20 keV (90% confidence level) and emission radius 75+ 3 , m for a distance of 0.5 kpc (implying 
a 0.5-10 keV luminosity ~ 6 x 10 29 erg s _1 ) is a viable interpretation, calling for a stream of parti- 
cles accelerated in the magnetosphere of PSR J0737-3039A and depositing their kinetic energy in the 
magnetic polar cap of PSR J0737-3039A or of the companion PSR J0737-3039B. A single power-law 
emission model implies a very steep photon index T = 4.2^ 2 'J, and a suspiciously high hydrogen column 
density whereas a photon index T = 2 does not provide an adequate description of the XMM-Newton 
and Chandra data. A two component model (a black body plus a power-law with r = 2) is statistically 
acceptable, but the additional power-law component is not required by the data. 
Subject headings: stars: neutron — pulsars: individual (PSR J0737-3039A, PSR J0737-3039B) - 
X-rays: stars — radiation mechanisms: thermal, non-thermal 



1. INTRODUCTION 

PSR J0737-3039A/B (Burgay et al. 2003; Lyne et al. 
2004) is the first double pulsar ever known. It comprises a 
22-ms pulsar (PSR J0737-3039A, hereafter A) and a 2.7- 
sec pulsar (PSR J0737-3039B, hereafter B) revolving in 
2.4 hr about the common center of mass along a mildly 
eccentric (e = 0.09) and highly inclined (~ 87°) orbit. 
These orbital parameters ensure a large magnitude for sev- 
eral post-keplerian parameters, making this binary an un- 
precedented test-bed for theories of gravity and relativistic 
physics (Lyne et al. 2004; Kramer et al. 2004). 

Further unique features of J0737-3039 double neutron 
star system are the occurrence of both an eclipse of A's ra- 
dio signal at superior conjunction (Lyne et al. 2004; Kaspi 
et al. 2004) and a modulation of flux and pulse shape 
of pulsar B along the orbit (Lyne et al. 2004; Manch- 
ester et al. 2004; Ramachandran et al. 2004). These 
phenomenologies are signatures of the interaction between 
the relativistic wind from A (releasing a spin-down power 
Ea = 5.8 x 10 33 erg s _1 , 3600 times larger than E B ) and 
B's magnetosphere and open the possibility of shedding 
light on still unanswered issues related with neutron star 
magnetosphere, pulsar winds and coherent radio emission 
(Lyutikov 2004; Arons et al. 2004; Jenet & Ransom 2004). 

The evidence for strong interactions between relativistic 
particles from one pulsar and the other's magnetic field im- 
mediately triggered interest for observations of the J0737- 
3039 system in the X-ray band. A short 10 ksec obser- 
vation with the Chandra ACIS-S instrument resulted in 
a clear detection of the binary (McLaughlin et al. 2004) 
with no significant evidence for variability along the or- 
bit on timescale longer than 15 min. The available 77 
photons could be represented satisfactorily (McLaughlin 
et al. 2004) by a power-law spectrum with photon in- 
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confidence error) resulting in a 0.2-10 keV luminosity of 
2.3±q;4 x 10 30 erg s" 1 at the distance of 0.5 kpc, inferred 
from the pulsars' dispersion measure and a model for the 
Galactic electron density (Cordes & Lazio 2002). The 
inferred X-ray luminosity and spectral index are consis- 
tent with an emission originating solely from the magne- 
tosphere of pulsar A, but emission from the shocked wind 
of pulsar A as it interacts with the interstellar medium 
(or with the magnetosphere of pulsar B) could be another 
viable hypothesis (Granot & Meszaros 2004). 

XMM-Newton targeted J0737-3039 system on April 
2004 collecting about three times more photons than the 
January 2004 Chandra pointing. We here report on the 
analysis of the XMM-Newton data and present the im- 
proved spectral fit made possible by joining the data set 
of the two observations. 

2. XMM-NEWTON OBSERVATION 

2.1. Data collection and reduction 

XMM-Newton observed PSR J0737-3039A/B on 10th 
April 2004 for 50 ks following an unsolicited target of op- 
portunity request. Data were taken in Full Frame (full 
image with 2.6 s timing resolution) and Small Window 
(~ 2' x 2' imaging with 0.3 s timing resolution) mode 
for MOS1 and MOS2, respectively, and with a medium 
filter. Data from the pn detector were taken in timing 
mode (1 dimensional imaging and 0.03 ms timing resolu- 
tion) with a medium filter too. We extracted the event 
files running emproc under SAS v6.0. Analysis of the light 
curve revealed that a sizable part of the observation was 
corrupted by soft photon flares. We extracted, from the 
entire MOS images, a total light curve with binning time 
of 50 seconds and created a time filter excluding intervals 
with rate larger than 5 count s _1 (e.g. Snowden et al. 
2004). The good time intervals amount to 20.5 ks and 
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25.4 ks for the two MOS cameras, respectively. In the 
central part of the MOS1 filtered image we can detect (us- 
ing XIMAGE v4.2) only one source. Its location is fully 
consistent with that of PSR J0737-3039A/B derived from 
timing observation (Lyne et al. 2004), with a count rate 
of (3.3 ± 0.5) x 10~ 3 couns s _1 . The source turns out to 
be too faint for a detection with the pn in timing mode 
and these data were not analysed any further. From the 
filtered data we extracted the MOS source spectrum (cir- 
cular region centered on source with r = 20" extraction 
radius) and the background spectrum (for the MOS1 us- 
ing an annular region centered on source with inner and 
outer radii ri nt = 40" and r ext = 80" respectively, and for 
the MOS2 using two r — 20" regions near the right corners 
of the small window image) . The background accounts for 
- 80% and - 65% of the total counts for the MOS1 and 
MOS2 respectively, which are 81 and 145. We generated 
response and ancillary files with the latest CCF library 
and fit together the two spectra. We have also analysed 
the data collected on 18th January 2004 with the ACIS- 
S instrument onboard Chandra finding similar results to 
McLaughlin et al. (2004). In the latter data, the source 
appears slightly extended (with a ~ 2a statistical signif- 
icance), but the small number of counts (about 80) does 
not allow to derive firm conclusions. 

2.2. Spectral fitting 

We tried several emission models to the joint MOS1, 
MOS2 and Chandra spectra rebinning the data to 20 and 
15 photons per channel for the two MOS and the ACIS-S 
respectively; hence a total of 16 spectral channels, ranging 
from 0.5 keV up to 10 keV for MOS1 and MOS2 (the pho- 
tons with energies less than 0.5 keV have been excluded 
owing to the uncertainties in the calibration of the two 
MOS cameras) and 0.3-10 keV for the ACIS-S data (cor- 
rected for the degradation in the ACIS-S quantum effi- 
ciency within XSPEC with the ACISABS component). The 
adopted spectral binning guarantees the proper use of the 
X 2 goodness-of-fit evaluator. The spectra have been mod- 
elled and fitted using the XSPEC package v. 11.3 with an 
absorption component modelled by TBABS. 

As reported in Table 1, all the explored single compo- 
nent models are acceptable from the point of view of the 
goodness-of-fit. However, their physical interpretation is 
less straightforward. We review these models in the fol- 
lowing. 

A black body model provides a good fit to the data 
with a xl e d = 1-6 (where xled * s the reduced x 2 )- Given 
the small number of degree of freedom (12) the null hy- 
pothesis probabily (n.h.p.) is 7% (Fig. 1). The black 
body equivalent temperature is fcTbb = O^Ot^Q 2 , keV (i.e. 

T = 2.5to'o5 x 10 6 K; here and in the following we use 
90% confidence interval for one interesting parameter, i.e. 
Ax 2 = 2.71). We can put only an upper limit on the col- 
umn density Nh < 0.4 x 10 21 cm~ 2 . This upper limit is 
consistent with the neutral hydrogen absorption measured 
in other sources with similar celestial coordinates and hav- 
ing distance ^ 1 kpc, as inferred using the ISM Column 
Density Search Tool 1 and from an estimate from the radio 
pulsars disperion measures of Njj = 1-5 x 10 20 cm -2 . The 
unabsorbed 0.5-10 keV flux is 1.9 x 10~ 14 erg s _1 cm~ 2 , 

1 http://archive.stsci.edu/euve/ism/ismform.html 



converting into a luminosity of 5.8 x 10 29 erg s _1 . For the 
reference distance of 0.5 kpc, the efficiency r] of conversion 
of rotational energy loss Ea into 0.5-10 keV luminosity 
is r\ ~ 10~ 4 and the equivalent black body radius is very 
small i?bb = 75lg m. Deriving an emitting radius of few 
kilometers requires much larger distances for the source, 
even when adopting neutron star atmosphere models. In 
particular, fixing the emitting radius to 10 km, the atmo- 
spheric emission model NSA of the XSPEC package places 
the source at 13.8 kpc, more than 20 times farther than 
the distance inferred according to the dispersion measure 
distance indicator. 

A thermal bremsstrahlung model is statistically accept- 
able too (xled = 1-0' n.h.p. 42%). The equivalent tem- 
perature is fcTbr = 0.56to;ig keV. The column density is 
Nh < 1.0 x 10 21 cm~ 2 The 0.5-10 keV unabsorbed flux is 
3.2 x 10 -14 erg s" 1 cm~ 2 . At 0.5 kpc this translates into 
a luminosity of 9.5 x 10 29 erg s . The emission measure 
for the adopted distance is 3.6tfl x 10 53 cm~ 3 . Following 
Grindlay et al. (2002) we can estimate the plasma density 
n <~ 10 26 i?~ 3 / 2 cm~ 3 for a spherical region of radius R. 
Assuming a mean particle density typical of the interstel- 
lar medium, i.e. n = 1 cm~ 3 , the radius of the emitting 
region turns out ~ 0.1 pc, much larger than the binary 
separation. Finally we note that a fit with a Raymond- 
Smith model provides an unacceptable result in the case 
of an emitting region of solar metalliticy (x 2 e d = 4-3)- The 
upper limit on the abundance is Z < 0.01 Zq. 

A power-law model (x 2 e d = 0-9, n.h.p. 54%) is charac- 
terized by a photon index T = 4.2~^ 2 '\. Given this steep 
power-law also the column density is very high Nh = 
1-8^0 7 x 10 21 cm~ 2 , comparable to the full galactic value 
(Nh = 4.7 x 10 21 cm -2 ) derived from neutral hydrogen 
measurements (Dickey & Lockman 1990). The unabsorbed 
0.5-10 keV flux amounts to 4.7 x 10~ 14 erg s _1 cm~ 2 . 
Placing PSR J0737-3039A/B at 0.5 kpc, this results in a 
luminosity of 1.4 x 10 30 erg s" 1 , corresponding to an effi- 
ciency T] ~ 2 x 10~ 4 . We note that aT = 2 model has a 
null hypothesis probability of 5 x 10~ 7 (see Fig. 2). 

A fit for a multi-components model keeping all the pa- 
rameters free is not meaningful, owing to the limited total 
number of available photons and taking in account the re- 
ported statistical significance of the fits for the single com- 
ponent models. Given its possible physical relevance (see 
below), we here only report on a two components model 
made by a black body plus a power-law with photon in- 
dex fixed to r = 2. Also in this case we obtained a good 
description of the data ( X 2 red = 1.2, n.h.p. 28%). The 
black body component parameters are similar to those of 
the single black body model (kT hh = 0.16 ± 0.04 keV, 
-Rbb = 90^28 m )- The percentage contribution to the total 
flux (2.7 x 10~ 14 erg s _1 cm~ 2 ) of the power-law compo- 
nent is ~ 52%. 

We also attempted a temporal analysis. We folded the 
MOS light curves to the orbital period without revealing 
any significant variation accross the orbit with an upper 
limit of 40% (3 a, see Fig. 3). The aforementioned failure 
in detecting the source in the pn detector prevented any 
search of modulation at the spin period of each the two 
pulsars. 
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3. DISCUSSION 

We have performed a joint spectral analysis of the 
XMM-Newton and Chandra observations of the J0737- 
3039 binary system, exploiting a three times larger amount 
of photons than the Chandra observation only (McLaugh- 
lin et al. 2004). 

Black body emission seems a viable model. Given the 
age of the system ( ^ 50 Myr since the second supernova 
explosion, Lorimer et al. 2004) the contribution in the 
0.5-10 keV band due to the cooling of the two neutron 
stars is negligible, but thermal emission could be powered 
by a flow of particles impinging onto the surface of one 
(or both the) neutron star(s). Pulsar A is likely to pro- 
vide the required energetic budget as the spin-down power 
of pulsar B is only 2 x 10 30 erg s _1 (Lyne et al. 2004). 
The simplest picture calls for back-flowing of charged par- 
ticles accelerated in the A's polar gap and depositing their 
kinetic energy in the A's magnetic polar cap (Cheng & 
Ruderman 1980; Arons 1981). The observed equivalent 
emitting radius is tiny, much smaller than the polar cap 
radius of pulsar A (r pc ~ 1 km), but partially filled po- 
lar cap thermal emission can be invoked, accounting for 
emitting radii as small as few tens of meters (Harding & 
Muslimov 2002). Alternatively, a much larger cap area 
results from modeling the data with a non-uniform tem- 
perature of the heated region (Zavlin et al. 2002). In this 
picture a sinusoidal X-ray emission at the pulsar A spin 
period (22.7 ms) is expected. 

In fact, a sizable number of millisecond pulsar (MSPs) 
with black body spectra has been already observed in the 
field and in globular clusters. In particular thermal con- 
tributions and cap radii as small as 30-40 m have been 
invoked in order to fit the X-ray spectra of the nearby 
pulsar PSR J0437-4715 (Zavlin et al. 2002) in the galac- 
tic field. We also note that the equivalent temperature of 
the emission from J0737-3039 system (0.2 keV) perfectly 
matches with the one of the pure black body spectrum 
of the MSPs observed in 47 Tuc's with Chandra (even if 
not based on a source by source spectral fits, Grindlay et 
al. 2002). The X-ray luminosities (and hence the emit- 
ting radii) are also comparable, once accounted for the 
intrinsic uncertainties in the distance of the J0737-3039 
system. This agreement is somewhat surprising given the 
much older ages of the MSPs in 47 Tuc (> few 10 9 yrs, 
Grindlay et al. 2002) with respect to the extimated time 
since formation of PSR J0737-3039A. 

Alternatively, it is conceivable that a fraction of the A's 
relativistic stream of particles is intercepted and channeled 
onto pulsar B polar cap(s) whose radius is 80 — 90 m (com- 
patible with the inferred Rbb at 95% confidence level). 
This hypothesis requires that ~ 3% of the spin-down lu- 
minosity (supposed isotropically released) of pulsar A pow- 
ers the observed X-rays. A modulation of the X-ray light 
curve at the B's spin period (or at a beating frequency 
between A's and B's rotational periods) is a predictable 
consequence of this model. A more detailed model along 
these lines has been presented by Zhang & Loeb (2004). 
Deeper exposures in small window with the pn instrument 



aboard XMM-Newton might reveal this modulation or the 
one at the A's spin period, strongly constraining the origin 
of X-ray in this binary. 

Bremsstrahlung emission requires an emitting region 
much larger than the binary separation and also larger 
(by a factor of ~ 100) than the radius of the bow shock 
between the relativistic wind of pulsar A and the ambient 
medium for a typical density of 1 atom cm -3 (and relative 
velocity of 200 km s _1 , Ransom et al. 2004). This trans- 
lates in an angular extension of the emitting region of ~ 1' 
for a distance of 0.5 kpc, which is excluded by both the 
Chandra and XMM-Newton observations. 

Absorbed single power-law emission models have been 
frequently used for interpreting the X-ray spectra of mil- 
lisecond pulsars (e.g. Becker & Trumper 1999; Becker & 
Aschenbach 2002; Nicastro et al. 2004). In this hypoth- 
esis the X-ray emission from J0737-3039 system would 
be ascribed to processes involving particles produced and 
accelerated in the magnetosphere of pulsar A. However, 
no MSP with an X-ray power-law spectrum as steep as 
r ^ 3 is known and a power-law model with a photon in- 
dex in the range 2 ^ T ^ 2.5 (the interval of values of T ob- 
served so far in the millisecond pulsar population, Becker 
& Trumper 1999) does not provide an adequate descrip- 
tion of the XMM-Newton and Chandra data. 

A simple power-law spectrum can naturally arise also 
from particles acceleration in a shock and this picture is 
particularly interesting in the case of the J0737-3039 sys- 
tem, where the relativistic wind from the recycled pulsar 
A could generate a bow shock at the interface with the 
insterstellar medium (provided the velocity of the binary 
through the medium is XL 200 km s _1 ) or when the wind 
interacts with the magnetosphere of pulsar B (Granot & 
Meszaros 2004). On the other hand, both these mech- 
anisms are believed to produce high-energy spectra with 
2 < r < 3 (Granot & Meszaros 2004) and a power-law 
index r £ 3 (at 90% confidence) appears difficult to acco- 
modate. 

A power-law spectrum with r = 2 can model satisfac- 
torily the observed spectrum of J0737-3039 system only 
if it is combined with a dominant thermal emission: in 
this case, at a distance of 0.5 kpc, the upper limit to the 
power-law contribution is at a level of ~ 4 x 10 29 erg s _1 . 
The latter luminosity would nicely match the predictions 
of Granot & Meszaros (2004) for the X-ray emission from a 
bow shock on the interstellar medium or surrounding pul- 
sar B magnetosphere. But also a simple magnetospheric 
emission from pulsar A 2 is compatible with that lumi- 
nosity. In fact, a two-component model (i.e. a thermal 
emission plus a non-thermal power-law contribution) has 
been adopted for improving the fit of the spectrum of PSR 
J0437-4715 (Zavlin et al. 2002), but only a long integra- 
tion with XMM-Newton can assess if it is required in the 
case of the J0737-3039 system. 

During the refereeing process we became aware of a sim- 
ilar analysis by Pcllizzoni et al. (2004). Even if they did 
not include Chandra data in the analysis they reached sim- 
ilar conclusions. 



2 It is calculated according to the correlation between the spin-down power E of the recycled pulsars in the galactic field and their observed 
X-ray luminosity L x , L x oc E@ with f} = 1.37±0.10, Possenti et al. (2002). The scatter around this correlation is about one order of magnitude. 
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Table 1 
Spectral fits to J0737-3039. 



Model Column density Ph. Index/ Xrod 
(10 21 cm' 2 ) Temperature (n.h.p.) 



Power-law 


i -°-0.7 


4 

*- Z -1.2 


0.9 (0.54) 


Power-law 


< 0.1 


2 (fixed) 


4.0 (1 x 10~ 6 ) 


Bremsstrahlung 


< 1.0 


0.56±g"i| keV 


1.0 (0.42) 


Black body 


< 0.4 


0.20 + ° : ° 2 keV 


1.6 (0.07) 


NSA 


< 0.5 


6.07« 


1.2 (0.24) 


Black+power 


< 0.9 


0.1618:8* kev 


1.2 (0.28) 



The degrees of freedom are 12 for all the models, 13 for the power-law with fixed r = 2 and 11 for the last model in the 
table. 
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Fig. 1. — Best fit of the pure black body spectrum. Squares indicate Chandra data, filled dots MOS1 data and open dots MOS2 data. 



7 




Fig. 2. — Countour plot of the power law fit, ellipses refers to 1, 2 and 3 a levels. A value of the photon index T = 2 is outside the 3cr 
contour level. Note that the column density axis is in logaritmic scale. 




Fig. 3. — MOS1 (upper) and MOS2 (lower) folded background subtracted light curves. Both curves are consistent with a constant emission. 



